Mean Field Theory of a Quantum Heisenberg Spin Glass 



Os 
ON 



C 
i 

C/3 



c3 



i 

c 

o 
o 



> ■ 

On ■ 

m : 

<N . 
On ■ 

O ■ 

ON : 
G\ . 

■ 



i 

c 

o 
o 



X 
S3 



Antoine Georges 1 , Olivier Parcollet 1,2 , and Subir Sachdev 3 
1 CNRS - Laboratoire de Physique Theorique, Ecole Normale Superieure, 24 Rue Lhomond 75005 Paris FRANCE 
2 Serin Physics Laboratory, Rutgers University, Piscataway, NJ 08854 USA 
3 Department of Physics, Yale University, P.O. Box 208120, New Haven, CT 06520-8120, USA 

(August 2, 1999) 

A full mean field solution of a quantum Heisenberg spin glass model is presented in a large- N limit. 
A spin glass transition is found for all values of the spin 5. The quantum critical regime associated 
with the quantum transition at 5 = 0, and the various regimes in the spin glass phase at high spin 
are analyzed. The specific heat is shown to vanish linearly with temperature. In the spin-glass 
phase, intriguing connections between the equilibrium properties of the quantum problem and the 
out-of-equilibrium dynamics of classical models are pointed out. 



The interplay between quantum effects and disorder in 
spin-glasses have been a subject of great recent interest 
[jlj. On the experimental side, the strength of quantum 
fluctuations can be continuously tuned by varying e.g. 
an applied transverse magnetic field (2j . Progress on the 
theoretical side, has followed two different routes. From 
the higher dimensional end, mean field solutions and ef- 
fective Landau theories have been obtained for quan- 
tum Ising and rotor spin glasses, with a special focus on 
the vicinity of the quantum-critical point where the glass 
transition temperature is driven to zero. In low dimen- 
sions f|,|],0] , it has been shown that the low T physics is 
controlled by rare events ( "Griffiths-McCoy" effects) at 
strong disorder fixed points. 

However, no established mean- field theory of the ex- 
perimentally important case of quantum Heisenberg spin- 
glasses, with full 577(2) symmetry, is yet available. Un- 
like the rotor/Ising models above, each site has non- 
trivial Berry phases which impose the spin commutation 
relations, and this is expected to place these models in a 
different universality class ||. Bray and Moore Q pio- 
neered the study of a model of Heisenberg spins on a fully 
connected ('Sherrington-Kirkpatrick') lattice of M sites. 
In this letter, we report a full solution of this model, both 
in the paramagnetic and the glassy phase, when the spin 
symmetry group is extended from SU(2) to SU(N) and 
the large-iV limit is taken. The Hamiltonian is 



(i) 



where the Jij are independent, quenched random vari- 
ables with distribution: P{Jij) oc e~ Ji i^ 2J \ In an imag- 
inary time path-integral formalism, the model is mapped 
onto a self-consistent single site problem with the action 
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S e n = S B + ^q j dTdT Q ( T - r')S a (r) ■ 5 6 (r') (2) 

with p = 1/ksT, and the retarded interaction Q ab (T — r') 
obeys the self-consistency condition : 



(r-r') = (l/N 2 )(s a (r)S b (r') 



Seff 



(3) 



Here, a, b = 1, • • • , n denote the replica indices (the limit 
n — > has to be taken later), and Sb is the Berry phase 
in the spin coherent state path integral. For N = 2 the 
problem remains of considerable difficulty even in this 
mean field limit. In Q, as well as in most subsequent 
work, the static approximation was used (see however 
ifjof), neglecting the t— dependence of Q afc (r); this may 
be appropriate in some regimes but prevents a study of 
the quantum equilibrium dynamics, in particular, in the 
quantum-critical regime. This imaginary time dynamics 
has however been studied recently in a Monte Carlo sim- 
ulation with spin S — 1/2 by Grempel and Rozenberg 
fLlfl , but their study was limited to the paramagnetic 
phase. In our large- A*" limit, the problem is exactly solv- 
able and, as explained below, this limit provides a good 
description of the physics of the N = 2 mean field model, 
as far as the latter is known. We find that in the param- 
agnetic phase, at low S (where the quantum fluctuations 
are the strongest), the quantum critical regime is a gap- 
less quantum paramagnet already studied in [p[p^] and 
radically different from the paramagnet obtained in the 
classical regime (at large 5), in which a local moment be- 
havior persists down to the glass transition. In the spin 
glass phase, various regimes are obtained as a function of 
temperature T. The thermodynamic properties and the 
dynamical response functions are analyzed below. Most 
notably, the low T specific heat is found to have a linear T 
dependence, a behavior commonly observed experimen- 
tally in spin-glasses but not often realized in mean-field 
classical models. Furthermore, the equilibrium dynamics 
of the quantum case reveals intriguing connections with 
some known features of the out-of equilibrium dynamics 
of classical glassy models, an observation already made 
in 1 13 1 in a different context. 

To handle the large- A limit, we use a Schwinger bo- 
son representation of the SU (A) spin operators : Sap = 
b a bp ~ SSa/3, corresponding to fully symmetric repre- 
sentations (one line of A^5 boxes in the language of 
Young tableaux) where the number of bosons is con- 
strained by J2 a b a b a — A5. In the SU(2) case, S co- 
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Quantum critical regime Quantum spin glass regime 

FIG. 1. Mean-field phase diagram and crossovers of the 
large- N quantum Heisenberg spin glass (the various regimes 
are discussed in the text). 

incides with the usual definition of spin. Fermionic rep- 
resentations can also be considered but they actually do 
not lead to a spin-glass phase at any temperature in 
the N = oo limit In the large- N limit, the self- 
consistent single-site problem reduces to a non-linear in- 
tegral equation for the replicated boson Green's function: 
G afc (r) = -£ Q (Tb°(T)bl b (0))/N§: 

{G- l ) ab {iv n ) = iu n 5 ab + X a S ab - ^ a \iv n ) (4) 
S Qb (r) = J 2 (G ab (r)) 2 G ab (-T) (5) 
G aa (r = CT) = -S (6) 

Here u n are the bosonic Matsubara frequencies, and G _1 
stands for the inverse in replica space. The (disorder- 
averaged) local spin correlation function is related to 

G ab (r) by xioc(r) = (^(0) • £,(r)) = G aa {r)G aa {-t) . 
The resulting phase diagram, obtained by both analyti- 
cal and numerical studies of these equations, is displayed 
in Fig. |], as a function of S and T / J. Spin-glass ordering 
is found at any value of S. The critical temperature in- 
creases as JS 2 at large S (see below) and vanishes in the 
limit £ — > 0, which corresponds to the quantum critical 
point in this model. Several crossovers are found within 
the spin-glass phase, which will be described later. 

We first describe the paramagnetic phase and the as- 
sociated crossovers. In this phase, the Green's func- 
tion is replica diagonal G ab (r) = G{r)5 a b and thus 
Eqs. (§-||) reduce to a single non-linear integral equa- 
tion. We emphasize that, as in any mean-field the- 
ory, paramagnetic solutions of the mean-field equations 
can be found even below the critical T where an insta- 
bility to ordering occurs. At high T, we have nearly 
free spins with an almost constant correlation function 
X\oc(t) — S(S + 1) and a Curie local susceptibility 
Xioc = Jo Xioc{r)dT ~ S(S + 1)/T. For large values of 
S, these solutions smoothly evolve, as T is reduced, into 
solutions which still behave locally as local moments, but 



with a Curie constant reduced by quantum fluctuations 
: xi oc ~ S 2 /T. This partial quenching occurs at a tem- 
perature of order JS 2 at large S, of the same order but 
smaller than the glass transition temperature. These so- 
lutions actually have unphysical low-T properties, such 
as a divergent internal energy U ~ — J 2 S 4 /2T and a neg- 
ative entropy (oc —J 2 S 4 /4T 2 ). These features are well- 
known in classical mean-field models and simply signal 
the tendency to spin-glass ordering. At smaller values 
of S (Fig.|l]) , a crossover to a different kind of paramag- 
netic solution is found below T ~ J, where we enter the 
quantum critical regime. In this gapless quantum para- 
magnet ('spin-liquid'), investigated previously in Hjjjl, 
the local response displays a scaling form for oj,T <C J, 
JX\ oc (oj) oc tanh(w/2T), and the local susceptibility di- 
verges only logarithmically Jxi oc oc ln( J/T). In contrast 
to the local-moment solutions, this paramagnet has fi- 
nite residual low-temperature entropy 0], so that the 
quenching of the entropy as T is decreased takes place 
much more gradually at low S when quantum fluctua- 
tions are strong, than at large S in the classical regime. 
It can be shown analytically E3] that these solutions of 
the mean-field equations exist down to T = only for 
very low values of S, smaller than S c — 0.05. For larger 
spins, a local- moment like solution is retrieved as T is 
lowered below a temperature of order JS 2 (again below 
the actual glass transition). However the spin-liquid solu- 
tions are the relevant ones in the quantum critical regime 
at finite temperature JS 2 < T < J for an extended range 
of spin values which extend up to S ~ 1. The detailed 
analysis of the coexistence between these two kinds of 
paramagnetic solutions at low S is rather intricate and 
will be presented elsewhere Jl3] . 

In the Quantum Monte Carlo results of |ll]] for the 
paramagnetic phase of the S — 1/2, SU(2) model, the 
same reduction of the Curie constant from S(S + 1) to 
S 2 was observed. Furthermore, the relaxation function 
x" {u))/oj evolves from a single peak of width JS centered 
at lo = to a three peak structure in the low-T local mo- 
ment regime. The central peak of weight S 2 corresponds 
to the residual local moment while two side peaks at an 
energy scale J 2 S 3 /T correspond to transverse relaxation 
|TTf . All these features are captured by our solution in the 
large- N limit, the only qualitative difference being that 
no thermal broadening of the central peak is found in this 
limit. Furthermore Jig] , numerical results not reported 
in fTl|| reveal that in a limited intermediate T range of 
the SU(2) 5=1/2 model, spin liquid solutions similar 
to those found here in the quantum critical regime are 
observed. Although a logarithmic regime is not directly 
visible in the T dependence of the local susceptibility be- 
cause of this limited range, quantum-criticality is directly 
apparent in a non-monotonic T dependence of the local 
spin correlation function Xlocij)- 

We now turn to the analysis of the spin-glass phase. 
We first note that the spin-glass transition is not signalled 
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by the divergence of the spin-glass susceptibility (which 
is actually of order 1/N) |lfj]. In the ordered phase, the 
boson Green's function can be parameterized as follows: 



G 



ab(r) = \ G(t) - (?) Sab ~ 9ab(l ~ S a b) 



(7) 



where g ab is a constant n x n matrix and g\ a constant, 
fixed so that G is regular at T = 0, i.e. G(t — ► oo) = 0. 
The usual spin-glass order parameter (l7) is q a b — 9a\- 
We have searched for replica-symmetry broken solutions 
by making a Parisi Ansatz for g a t, and found that single- 
step replica symmetry breaking always applies []l0|| . The 
Parisi function g(x) associated with g a {> is thus piecewise 
constant: g(x) = for x < x c , g{x) — g(l) = ^fq~EA = g 
for a; > x Cl where qEA is the Edwards- Anderson order 
parameter; this also implies that g — g. For the following 
discussion, it is convenient to define the parameter = 
—JG(iv — 0)/g. Using standard inversion formulas for 
a Parisi matrix |l8|, the full set of mean- field equations 
read: 

(G(ii/ n )) ~* = iv n - Jg/Q - (£(»!/„) - S(0)) (8) 
E(r) = J 2 (g 2 (t)G(-t) - 2 5 G(r)G(-r) 



G(r = (T) 

j3x c 



-gG 2 (T)+2g 2 G(T)+g 2 G(-T)) (9) 



g — S 

(i/e-e)/j 9 2 



(10) 

(ii) 



However, these equations do not determine O (or equiv- 
alently the breakpoint x c ) as also happens in classical 
spin-glass model with a single-step of replica symme- 
try breaking: there is a continuous family of solutions 
parametrized by 0, which has to be determined by inde- 
pendent considerations. Two possibilities have appeared 
in previous work: (i) Determine by minimizing the 
free energy, as a function of 0, or (ii) impose a vanishing 
lowest eigenvalue of the fluctuation matrix in the replica 
space (the "replicon" mode). Criterion (?) is certainly 
the natural one from the point of view of equilibrium 
thermodynamics. However, studies of out-of-equilibrium 
dynamics of classical spin-glasses have revealed jl9| that 
these lowest free-energy solutions can never be reached 
and that the system "freezes" at a dynamical tempera- 
ture T° , given precisely by the onset of solutions satisfy- 
ing the "replicon" criterion (ii). In our quantum problem, 
both choices give sensible solutions, but with entirely dif- 
ferent spectra of equilibrium dynamical fluctuations: (i) 
leads to a gap in XiocC 1 -^)' wnuc (**) is found to be the 
unique choice leading to a gapless spectrum. A similar 
observation was made in the work of Giamarchi and Le 
Doussal in their study of a one-dimensional quantum 
model with disorder. In the present context, it seems nat- 
ural to expect local gapless modes in the ordered phase 
of a quantum spin-glass with continuous spin symmetry, 
and these various considerations lead us to adopt (ii). 



Diagonalizing the fluctuation matrix in replica space, we 
find the lowest eigenvalue ei = 3/?J 2 g 2 (l — 30 2 ) so that 
the replicon criterion leads to = l/v3> independent 
of T; the same value also appears by independently im- 
posing that G has a gapless spectral weight. In contrast, 
criterion (i) leads to 21n0+l/(40 2 ) + l/2-30 2 /4 = 0, or 
~ 0.44 . . ., and a gapped solution. We also note that 
the previous computation shows that the replica sym- 
metric solution = 1 is unstable in the spin glass phase. 
Moreover, it can be shown that it leads to unphysical 
negative spectral weight at large-S*. Hence, a correct de- 
scription of the low-energy excitations of the quantum 
model requires replica symmetry breaking at any finite 
T in the spin-glass phase, although the replica symmetry 
is restored at T = where x c = (from (|TT|)). 

Once is determined, a full numerical solution of 
the above equations can be performed. In particular, 
the "equilibrium" spin glass temperature T s e g 9 obtained 
from criterion (i) is lower than the "dynamical" tran- 
sition temperature T s c g obtained from criterion (ii) (see 
Fig. |IJ): this is not obvious a priori, but is certainly re- 
quired in our interpretation. Further analytical insight 
can be obtained in the limit of large S. This limit can 
actually be taken in two distinct ways, revealing two 
crossovers within the spin-glass phase displayed in Figjl]. 
If we take S — > oo while keeping T/JS 2 fixed (i.e stay- 
ing close to the critical temperature), all non-zero Mat- 
subara frequencies can be neglected (the static approxi- 
mation is accurate). In this limit, we find in particular 
T s c g - 2JS 2 /3 3 / 2 . Alternatively, keeping f = T/JS and 
Q = oj/JS fixed, we access the "semi-classical" regime 
of the spin-glass phase. In this limit, the Green's func- 
tion obeys a scaling form G(uj,T) = f(Cj)/(JS), where / 
turns out to be independent of T and satisfies : 



/(Q)- 1 = u, - 1/0 - 30 - f(Q) - f*(-u>) 



(12) 



Eliminating /*(—£)) leads to a quartic equation for f(ui) 
on which all the above properties can be checked more 
explicitly. A plot of the (gapless) relaxation function 
in the spin-glass phase x"(w)/w obtained from ( |l~2| ) is 
displayed in Fig.||. 

Finally, we briefly describe the thermodynamic prop- 
erties, focusing on the T dependence of the specific 
heat. Numerical results for this quantity for interme- 
diate spin are displayed on Fig ||. They have been 
obtained from the T-derivative of the internal energy 
U = -J 2 /2 f£ G ab (T) 2 G a b(-T) 2 dT, where G is a nu- 



merical solution of Eqs.(|8Hll|). Furthermore, a large-.!?, 
low— T expansion of U(T) can be done analytically and 
leads to 0: U(T) = U(0) + aST 4 + bf 2 + ... where a 
and b are positive numerical coefficients. Hence in the 
quantum regime defined by T < J\/~S (see Fig. 0), the 
specific heat depends linearly on temperature. Moreover, 
this behavior actually holds numerically for intermediate 
values of the spin as displayed in Fig |^. 
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Note added: It has been recently proven by one of us 
pl[ that the behaviour Jx" oc (w) oc const, found above 
in the quantum critical regime also holds for the SU (2) 
case. 



FIG. 2. Relaxation function x"^)/^ m the large-5* limit, 
obtained from (p^). 




FIG. 3. Specific heat C(T) and internal energy U(T) (in- 
set) vs. temperature T, from a numerical solution of Eqs. 
(10) for 5 = 5. 



Despite being formulated over two decades ago j9|, a 
complete understanding of the quantum Heisenberg spin 
glass at the mean-field level has proven elusive. Here we 
have obtained a complete solution in a large N limit, and 
presented evidence that global aspects of the phase dia- 
gram pertain also to the physical SU(N — 2) case. We 
described crossovers in the vicinity of a quantum critical 
point accessed by varying the spin S, but we can expect 
that some features and intermediate temperature regimes 
will survive when it is accessed by varying other param- 
eters in the Hamiltonian, including doping with metal- 
lic carriers as in Kondo lattice models |2(|[l2|]. We have 
also described the T — > thermodynamics and spectral 
functions within the spin-glass phase, which is something 
not previously analyzed in any mean-field quantum spin 
glass model: we found a specific heat linear in tempera- 
ture, and a dynamical susceptibility x"(a;)/u> — > const as 

LU -> 0. 
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